Mem. S.A.It. Vol. 0, 1 

©SAIt 2008 



Memorie deiia 




Uegli S9° 



Granulation across the HR diagram 

I. Ramirez^ C. Allende Prieto^, D. L. Lambert^, L. Koesterke^-^, and M. Asplund 

' Max Planck Institute for Astrophysics, Postfacli 1317, 85741 Garching, Germany 

^ Mullard Space Science Laboratory, University College London, UK 

^ McDonald Observatory and Dept. of Astronomy, University of Texas at Austin, USA 

Texas Advanced Computing Center, University of Texas, USA 

e-mail: ivan@mpa-garching.mpg.de 



Abstract. We have obtained ultra-high quality spectra (R = 180,000; S /N > 300) with 
unprecedented wavelength coverage (4400 to 7400 A) for a number of stars covering most 
of the HR diagram in order to test the predictions of models of stellar surface convection. 
Line bisectors and core wavelength shifts are both measured and modeled, allowing us to 
validate and/or reveal the limitations of state-of-the-art hydrodynamic model atmospheres 
of different stellar parameters. We show the status of our project and preliminary results. 
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1. Motivation 

Convection is one of the least understood phe- 
nomena in stellar astrophysics. This is, in part, 
due to the lack of high quality data that can 
be used to test theoretical models. Such tests 
can allow us to better understand the physics 
involved in the convection phenomenon and 
may help us to reveal the limitations of state- 
of-the-art models. In addition to this very fun- 
damental motivation, we note that an accu- 
rate determination of abundances and funda- 
mental parameters relies on models providing 
an adequate description of the outer layers of 
stars. Unrealistic models may lead to uncertain 
abundances and these errors can easily propa- 
gate to our interpretations of Galactic chemi- 
cal evolution as well as stellar and primordial 
nucleosynthesis. The impact of granulation on 
stellar abundances has been shown to be non- 
negligible in many cases (e.g., Asplund i2005i) . 
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2. Sample and strategy 

Our sample consists of 15 stars cover- 
ing a substantial part of the HR dia- 
gram and includes one metal-poor giant. 
Observations are being completed with the 
Robert G. TuII coude spectrograph on the 
2.7 m Harlan J. Smith Telescop e at McDonald 
Observatory jTuII et af] 1 19951) . The average 
spectral resolution of our data is 180,000 and 
signal-to-noise ratios range from about 300 to 
500. The spectra cover (completely in the ma- 
jority of cases) most of the visible range (4400 
to 7400 A) being thus the most complete data 
set of this quality collected to date. 

Line bisectors and core wavelength shifts 
are measure d for u nblended Fe i features listed 
in the iNave et all (I'l994l) catalog. Line bisec- 
tors of features of similar line-depth are aver- 
aged to reduce the noise given that the shape of 
granulation bisectors is determined mainly by 
the line strength. These average hne bisectors 
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Fig. 1. Average bisectors of Fe i lines sorted in four groups of similar line-strength for 4 of our sample 
dwarf stars and the subgiant star Procyon, covering roughly the temperature range from 4800 to 6600 K. 
The cores of these bisectors have bee n set to be at z ero velocity in order to reduce the noise introduced by 
uncertain laboratory wavelengths (cf. lRamfrez et al ]|2b08ll200Si) . 



are then compared to their theoretical counter- 
parts computed from 3D hydrodynamic mod- 
els. If the comparison shows reasonably good 
agreement, the model is considered valid and 
can be reliably used in the derivation of 3D cor- 
rections to stellar parameters and abundances. 

We compute a few lines that have very ac- 
curate laboratory data available using both ID 
and validated 3D models to estimate 3D cor- 
rections to chemical abundances derived with 
the standard ID approach. We also compute 
the full spectrum using extensive atomic and 
molecular line databases to look for features 
that are potentially sensitive to 3D effects as 
well as to explore the impact of surface in- 
homogeneities on absolute fluxes (for exam- 
ple to study their impact on synthetic colors). 
The necessary tools to compute spectra from 
3D hydrodynamic models have already been 
developed (e.g., Koesterke et al. 2008). 

As of this writing, we have successfully 
applied the procedure described above to a 
samp le of K-dwarf stars (Ramirez et al. 2008, 
120091) . We are currently working on an ex- 
tension to stars of different stellar parameters. 
Most of the needed observational material is al- 
ready at hand. Below, observational results for 
a few dwarf stars are briefly summarized. 



3. Results: the main-sequence 

Preliminary results for 5 of our sample stars 
are shown in Figs.[T]and|2] These stars roughly 
cover the main-sequence from mid-F to early- 
K spectral type (note however that Procyon is 
a subgiant). 

The trends observed are in good agree- 
ment with tho s e found by o ther authors (e.g., 
iDravinsI Il987l: iGrayl l2005l) . The line bisec- 
tor span increases with efi'ective temperature; 
those for the strongest lines go from less than 
100 ms ' in the K-dwarf HIP 86400 to about 
0.6 kms ' in Procyon (Fig. [T]i. Weaker lines 
show smaller bisector span. The bisectors have 
the typical C-shape that is associated with the 
granulation phenomenon. 

Core wavelength shifts behave similarly. 
Convective blueshifts are larger (i.e., wave- 
length shifts are more negative) for the weaker 
lines of warmer stars (Fig. 121). This is due to 
the fact that such lines form in deeper atmo- 
spheric layers where the granulation contrast is 
larger and because convection is less efficient 
in warmer stellar atmospheres, which means 
that larger velocities are necessary to transport 
the convective energy. 
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Fig. 2. Left panels: Core wavelength shifts as a function of equivalent width (EW) of Fe I lines in the 
spectra of four of our sample dwarf stars and the subgiant star Procyon, covering roughly the temperature 
range from 4800 to 6600 K. Right panels: As in the left panels but with a logarithmic ordinate axis scale. 
Dots correspond to individual lines; filled circles with error bars represent bin-averaged values. The range of 
convective blueshifts increases from about 200 to 1000 m s"' from the K-dwarf HIP 86400 to the F-subgiant 
Procyon. The blueshifts decrease for stronger lines but reach a saturation value (a plateau best seen in the 
left panels, forced to be at zero velocity) that begins at systematically higher EW for warmer stars. Note 
also that a plateau seems to be present at the low EW end (best seen in the right panels). 
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Along the main-sequence, the effective 
temperature seems to be the only parameter 
determining the shapes of line bisectors and 
the core wavelength sh ift vs. line strength rela- 
tions but note that Ramirez et al.l t2008, 20091), 
who studied these effects in K-dwarfs, have 
found that stars showing high levels of chro- 
mospheric activity presented anomalies. The 
very active star e Eri, for example, shows al- 
most vertical line bisectors and a very noisy 
core wavelength shift vs. EW relation com- 
pared to other K-dwarfs even though its spec- 
trum is of the highest quality. 

In Fig. [3] we compare the core wave- 
length shift vs. EW relati on obtained us ing 
the solar spectrum by Kuru cz et al.l (1 19841) as 
well as our skylight observations with that 
corresponding to the solar-twin star 18 Sco 
jPorto de Mello & da Silval Il997h . The latter 
has fundamental parameters essentially iden- 
tical to those of the Sun and thus one would 
expect it to show identical granulation signa- 
tures. However, it is clear that the convective 
blueshifts for the weakest lines in the 18 Sco 
spectrum are at least 50 ms"' too small (less 
negative) compared to those for the Sun. 

It is interesting to note that the data for 
18 Sco were obtained at an epoch when the 
star was experiencing a maximum of activity 
(May of 2007), s ignificantly stronger than the 
Sun's (iHall et al.ll2007i) . Whether some activity 
related effect is responsible for this difference 
needs to be confirmed by follow-up observa- 
tions of the star during a minimum of chro- 
mospheric emission along with nearly simul- 
taneous (and homogeneous; i.e., same spectro- 
graph) observations of a very bright asteroid or 
other sunlight-reflecting object. 

4. Summary 

We have been able to detect, quantify, and 
successfully model the granulation signatures 
present in K-dwarf spectra. Moreover, we have 
used our validated simulation to explore the 
impact of 3D effects on stellar abundances and 
fundamental parameters. Now we plan to ex- 
tend this work to stars of different stellar pa- 
rameters covering most of the HR diagram, in- 
cluding a metal-poor giant. 
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Fig. 3. Bin-averaged core wavelength shift vs. 
equivalent width relation (cf. Fig. IS for the solar 
spectrum by Kurucz et al. (1984), degraded to a 
spectral resolution similar to that obtained in our 
work (red dotted line) as well as our skylight (green 
dashed line) and 18 Sco (black solid line) observa- 
tions. 

Preliminary results for dwarf stars confirm 
previous findings that the granulation effects 
are more important for warmer stars. We find 
that the effective temperature is possibly not 
the only parameter determining the strength of 
these signatures as the comparison of obser- 
vational results for active and inactive K-type 
dwarf stars of essentially the same temperature 
revealed important differences. Furthermore, 
the granulation signatures of the solar-twin star 
18 Sco at its activity peak differ significantly 
from those observed in the Sun. 
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